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Abstract 

By employing the higher (N>5) dimensional version of the Wu-Yang Ansatz we obtain magnet- 
ically charged new black hole solutions in the Einstein- Yang-Mills-Lovelock (EYML) theory with 
second (02) and third (as)order parameters. These parameters, where 02 is also known as the 
Gauss-Bonnet parameter, modify the horizons (and the resulting thermodynamical properties) of 
the black holes. It is shown also that asymptotically (r — > 00), these parameters contribute to 
an effective cosmological constant -without cosmological constant- so that the solution behaves 
de-Sitter (Anti de-Sitter) like. 
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I. INTRODUCTION 

As the requirements of string theory/brane world cosmology, higher dimensional (N > 4) 
space times have been extensively investigated during the recent decades. Extensions of 
N = 4 Einstein's gravity has already gained enough momentum from different perspectives. 
General relativity admits local black hole solutions as well as global cosmological solutions 
such as de Sitter (dS) and anti de-Sitter (AdS), which are important from the field the- 
ory (i.e., AdS/CFT) correspondence point of view. Once a pure gravity solution (with or 
without a cosmological constant) is found the next (routine) step has been to search for the 
corresponding Einstein-Maxwell (EM) solution with the inclusion of electromagnetic fields. 
As a result EM dS (AdS) space times have all been obtained and investigated to great ex- 
tend. One more extension from a different viewpoint, which is fashionable nowadays, is to 
consider extra terms in the Einstein-Hillbert action such as the ones considered by Lovelock 
in higher dimensions to obtain solutions, both black holes and cosmological [1] . The added 
extra terms in the action have the advantage, as they should, that they don't give rise to 
higher order field equations. It is known that in N = 4, the second order (also known as 
Gauss-Bonnet) Lovelock Lagrangian becomes trivial unless coupled with non-trivial sources 
such as non-minimal scalar fields [9]. In higher dimensions (N > 5), however coupling with 
electromagnetic fields proved fruitful and gave rise to interesting black hole/cosmo logical 
solutions[2]. In this regard, Brihaye, et al, have worked on particle-like solutions of EYM 
fields coupled with GB gravity in N-dimensional spherically symmetric spacetime [5]. To 
have a non-trivial theory with the higher order Lovelock Lagrangian with third order pa- 
rameter on the other hand, we need the dimensionality of our space time to be N > 7. In 
this paper we shall follow similar steps, to extend the results of electromagnetic fields to the 
Yang-Mills (YM) fields with gauge group SO(N — 1). As expected, going from Maxwell to 
YM constitutes a highly non-trivial process originating from the inherent non-linearity of 
the latter. To cope with this difficulty we employ a particular YM Ansatz solution which was 
familiar for a long time to the high energy physics community. This is the Wu-Yang Ansatz, 
which was originally introduced in N = 4 field theory [3, 4]. Recently we have generalized 
this Ansatz to N = 5, in the Einstein-Gauss-Bonnet (EGB) theory and obtained a new 
Einstein- Yang-Mills-Gauss-Bonnet (EYMGB) black hole [5]. By a similar line of thought 
we wish to extend those results further to iV > 5 and also within the context (for N > 7) 
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of third order Lovelock gravity. Our results show that both the second (a 2 ) and third (a 3 ) 
order parameters modify the EYM black holes as well as their formation significantly. For 
instance in the g tt term the gauge charge term comes with the opposite sign and fixed power 
of \ , which is unprecedented in the realm of EM black holes of higher dimensions. This 
makes construction of black hole types from pure YM charge (with negligible mass, for ex- 
ample) possible and enriches our list of black holes with new properties. What follows for 
N > 7, for technical reasons, we assume a relation between a 2 and a 3 which are completely 
free otherwise. Not only black holes but the asymptotical behaviors and properties of our 
space times are determined by these parameters as well. It is not difficult to anticipate 
that by the same token the same problem can further be generalized to cover the forth (014), 
fifth (0:5), etc. order terms in the action to be superimposed to the Einstein-Hillbert (EH) 
Lagrangian. By studying the relative weight of contribution from higher Lovelock terms it is 
not difficult to anticipate that the EH and GB terms dominate over the higher order, much 
more tedious terms. For this reason we restrict ourselves in this paper to maximum third 
order (a 3 ) terms in the Lagrangian. 



II. ACTION AND FIELD EQUATIONS 

The action which describes the third order Lovelock gravity coupled with Yang-Mills field 
without a cosmological constant in N dimension reads [1] 

la = \ [ dx N ^Tg [C EH + a 2 C GB + a 3 C (3) - tr (F$F^)] , (1) 

(JV-l)(JV-2)/2 

where tr{.) = Yl (•) > £eh — R is the Einstein-Hilbert Lagrangian, Cgb = 

a=l 

R^ujsR^ 16 - ^R^BT + R 2 is the Gauss-Bonnet (GB) Lagrangian, and 

f O TDUVCTK p TDpT I Q TDCTK T3pT 

M3) - m tiunprri ^ + Sit cp K UT K MK 

I O/ p/l" t " : D DP I QDDF'" : D 

+ Z^JX ti UKUp ti ^ + otiti riannv 

+ 2AR^ K R^R KV + \m^R vu R\ (2) 
— Y1RR^ P R^y + R^ 1 

is the third order Lovelock Lagrangian. Here R, R^5 and R^ u are the Ricci Scalar, Riemann 
and Ricci tensors respectively, while the gauge fields F$ are 

p(a) _ a A{a) _ o A(a) , _}_M a ) /l( fe )/l( c ) (<1\ 



where C^j^ are the structure constants of ( JV ~ 1 K JV ~ 2 ) -parameter Lie group G, a is a coupling 
constant, are the gauge potentials, and a 2 and a 3 are GB and third order Lovelock 
coefficients. Variation of the action with respect to the space-time metric g pv yields the 
Einstein- Yang-Mills-Gauss-Bonnet-Lovelock (EYMGBL) equations 

G% + «2Ggf + a 3 G$ = T^, (4) 
where the stress-energy tensor is 



= tr 



(5) 



G pu is the Einstein tensor, while C^f and G^ 3 J are given explicitly as[2] 

G>if = 2 {—Rii,<j KT R KT(J v — 2R lxpi/cr R pa — 2R fia R J u + RR pu ) — -Cgb9pv, (6) 

Gjg = -3 (4it^™ R aK \pR^ l/Tfl - 8R Tp Xa R™ Tfl R x upK (7) 

I Q P TCTK D F?^/ 9 _ DTpCTK p p I O pr P CTK PP 

- r^- f \< J^anXp^ T fi ^anrp^vp ~T~ Oil yap 1 ^ tu k 

_|_Q pO" t? T P t? K I /ID TCTK p pp zip TCTK p pp 

+AR TpCK RcKrpRvp + 2RR u RTp R TpKll + 8R T upp R p a R p T 

+4R Tp R pr R vp — 8R T l/ R Tp R p ll + 4RR I/p R p fl 
-R 2 R up ) - -£(3)57^. 

Variation of the action with respect to the gauge potentials yields the Yang-Mills 
equations 

+ icW (c) AWF^ = 0, (8) 

while the integrability conditions are 

♦jrW^ + ±C?g> c) 4» * F (c) " v = 0, (9) 

in which * means duality [6]. 



III. WU-YANG ANSATZ IN N > 5 DIMENSIONS 

The N— dimensional line element is chosen as 



ds 2 = -f(r) <& + Jfi + r " d & N - 2 , (10) 



in which the S N 2 line element will be expressed in the standard spherical form 

AT-3i-l 

dn%_ 2 = d9l + E II sin 2 6j d9l 0<e 1 <n,0<e i <2n. (11) 

i=2 j=l 

We use the Wu-Yang Ansatz [7] in N— dimensional case as 

Q N -* 

= — (xidxj — Xjdxi) , Q = charge, r 2 = y^x 2 , (12) 
^ i=i 

2<j + l<i<N-l, and 1 < a < (N - 1) (N - 2)/2, 

where we imply (to have a systematic process) that the super indices a is chosen according 
to the values of % and j in order. 

The YM field 2-forms are defined as follow 

F (a) = dA (a) + 1 Ma) Mb) A ^(c) _ ( ^ 

We note that our notation follows the standard exterior differential forms, namely d 
stands for the exterior derivative while A stands for the wedge product [7]. The integrability 
conditions 

dF^ + ^.C$ {c) A^ A = 0, (14) 
are easily satisfied by using (12). The YM equations 

d * F {a) + ^C$ ic) A<V A *F (c) = 0, (15) 

are also satisfied. The energy-momentum tensor (5), becomes after 



tr 



F (a) p(a)\a 



(N — 3) (iV-2) Q 2 

r 4 



(16) 



with the non-zero components 



H= (JV - 3) ^- 2)<y diag [-1,-1, «,«,..,«], and «=-^|. (17, 

The EYMGBL equations (4) reduce to the general equation 

(r 5 - 2a 2 r 3 (/ (r) - 1) + 35 3 r (/ (r) - l) 2 ) /' (r) + 
(n - 1) r 4 (/ (r) - 1) - (n - 3) a 2 r 2 (/ (r) - l) 2 + 

(n - 5) a 3 (/ (r) - l) 3 + (n - 1) r 2 Q 2 = 0, (18) 
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in which a prime denotes derivative with respect to r, n = N — 2, a 2 = (n — 1) (n — 2) a 2 
and o 3 = (n — 1) (n — 2) (n — 3) (ra — 4) o 3 . This equation is valid for N > 4 (i.e. n > 2), 
but for TV = 4 (i.e. n = 2) we get 

r 3 /'(r)+r 2 (/(r)-l) + Q 2 = 0, (19) 

which clearly is 0:2,3 independent and therefore it will be the Einstein- Yang-Mills equation 
admitting the well-known Reissner-Nordstrom form 

For N = 5 (i.e. n = 3) the Eq. (18) has already been considered in [5]. We note that in 
Eq. (18), r, a 2 and a 3 can all be scaled by Q (i.e., r — > |Q| r, o 2 — * Q 2 o 2 and o 3 — > Q 4 a3, 
as long as Q 7^ 0) so that we set in the sequel, without loss of generality, Q — 1. 



IV. THE EYMGB CASE, a 2 + 0, a 3 = 

Eq. (18) with a 3 = takes the form 

(r 3 -2a 2 r(f(r)-l))f'(r)+ (21) 
(n - 1) r 2 (/ (r) - 1) - (n - 3) a 2 (/ (r) - l) 2 + (n - 1) = 0, 

which may be called as Einstein-Gauss-Bonnet- Yang-Mills (EGBYM) equation. This equa- 
tion admits a general solution in any arbitrary dimensions N as follows 

where m is the usual integration constant to be identified as mass. 



A. The EYMGB solution in 6— dimensions 



In this section we shall explore some physical aspects of the solution (22) in 6-dimensions. 
This is interesting for the reason that, N = 5 and N = 6 are the only dimensions which will 
not be effected by the non-zero third order Lovelock gravity. For N = 6 (n — 4), the metric 
function f(r) in Eq. (22) takes the form 



/± W = 1 + ^- 1 1 ± y 1 + + —pr I . (23) 



in which a 2 (— Gct 2 ) and ± refer to the two different branches of the solution. Asymptotic 
behaviors of f± (r) can be shown to be as 

lim f + (r) — > 1 + — , and lim /_ (r) — > 1, 

which imply that, the positive branch is Asymptotically-de Sitter (A-dS) with positive a 2 
and Asymptotically-Anti de Sitter (A-AdS) with negative a 2 . It is seen obviously that the 
negative branch is an Asymptotically Flat (A-F) space. One can also show that 

lim / + (r) — > +00, and lim /_ (r) — > —00, 

r— >0+ r^0+ 

which clearly, shows that, /+ (r) is an A-dS solution while /_ (r) represents an A-F black 
hole solution. In the sequel we shall consider a 2 > with the negative branch of the solution 
(i.e., the A-F black hole solution). One can easily show that, this solution admits a single 
horizon (i.e., event horizon) given by 

"+ = (t + W 7 ^) " + (t - W 7 ^) " • < 24 > 

which is real and positive for any values of m and a 2 - In Fig. (1), (i.e., the dashed curves) we 
plot the radius of event horizon r + , in terms of ^(i.e., Gauss-Bonnet parameter), for some 
fixed values for m. This figure displays the contribution of the Gauss-Bonnet parameter to 
the possible radius of the event horizon of the black hole. By looking at the Fig. (1), one 

may comment that for any value of m, lim r + — > 0. We notice further that, with a 2 = 0, 

02^00 

we get the radius of event horizon for the six dimensional EYM black hole solution which 
was given in the Ref. [5]. That is, for very large a 2 , the event horizon coincides with the 
central singularity. As a particular choice, we consider a 2 — \ which implies that 

r+ = m 1/3 . (25) 



The surface gravity, k defined by 



^ = -\g tt g ij 9tt,i 9u tj , (26) 



takes the value 

' 1 



K 



/(r + ) 



m l/3 



2m 2 /3 + 6' 

The associated Hawking temperature depending on mass m and a 2 — \ becomes 



(27) 



T K 3 

Th ~ 27 _ 4^m2/3 + 6 ' ( 8) 



in the choice of units c — G — h—k — 1. 
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B. The EYMGB solution in 7-dimensions 



In this section we represent some physical aspects of the solution (22) in 7-dimensions. In 
7-dimensions, both second and third order Lovelock terms contribute but still we set a 3 = 
in order to identity the contribution of a 3 . For N = 7 (n = 5), the metric function /(r) in 
Eq. (22) takes the form 

n / n r 2 I / Ad^rn 8a 2 \ , s 

f *v =1+ 2*( 1± r + — + -?r)' (29) 

in which a 2 (= VZot-i) and ± refers to the two individual branches of the solutions. In Fig. 
(2) we plot /_ (r) which goes to asymptotically flat value for r — > oo. Asymptotic behaviors 
of f± (r) can by written as 

lim / + (r) — > 1 + — , and lim /_ (r) — > 1, 

r^oo CK2 r-+oo 

which imply that, the positive branch is Asymptotically-de Sitter (A-dS) for positive a 2 and 
Asymptotically- Anti de Sitter (A-AdS) for negative a 2 , while the negative branch leads to 
an Asymptotic Flat (A-F) space. Also one can show that 

lim / + (r) — > +oo, and lim /_ (r) — > — oo, 

which clearly, manifests that, /+ (r) is an A-dS solution while /_ (r) represents an A-F black 
hole solution. Hence in the sequel we just consider a 2 > and the negative branch of the 
solution (i.e., the A-F black hole solution). One can easily show that, this solution admits 
only an (event horizon) which can be written as 

r+ = \jl - 6a 2 + \J(l - 6a 2 ) 2 + 4^, (30) 

which implies that r + is real and positive for any values of m and a 2 . In the Fig. (1) we 
plot the radius of event horizon r + , in terms of a 2 (i.e., the solid curves), with some fixed 
values of m. This figure displays the contribution of the Gauss-Bonnet parameter in place 
of possible radius of the event horizon of the EYM black hole. We notice that, with a 2 = 0, 
we recover the radius of event horizon for the 7— dimensional EYM black hole solution which 
is given in the Ref. [7]. 

As a particular choice, we consider a 2 = 1/6 which implies that 

r + = (4m) 1/4 , (31) 
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and the surface gravity, (26) in this case has the value 



K 



(32) 



m l/4 + 4- 

With the associated Hawking temperature given by 

K \ Tfl} ^ 4 

Th = 2^ = nm^Tl- (33) 
which is comparable with the 6-dimensional case (28). 

V. THE EYML CASE WITH, a 2 = 0, a 3 + 

In this section we just consider the effect of a 3 on the solution of the field equation. Eq. 
(18) with a 2 = takes the form 

(r 5 + 3a 3 r (f (r) - l) 2 ) /' (r) + z/r 4 (/ (r) - 1) + (34) 
(z/-4)a 3 (/(r)-l) 3 + z/r 2 = 0, 

in which v = N — 3 (= n — 1 therefore). The proper general solution of this equation is 
given by 

/(r) = 1+ (35) 

6a 3 (i/-2)r 3 fts 

in which 



Q = -108A + 12W — [275 3 A 2 + 4 (i/ - 2) 2 r 2 ^+ 4 l , 

V «3 

and 

A = vr v ~ 2 + (y- 2)m. 
One can easily show that in the limit <S 3 — > 0, we obtain 



m v 
&3^o J v ' " r v {v — 2) r 

which is the EYM solution. 



lim/(r) = l--- ,., (36) 



By setting v = 4 in Eq. (39) for N = 7 one obtains 



Qs 4r 4 
12a 3 



in which 



fl — -108 (4r 2 + 2m) + 12 [27a 3 (4r 2 + 2m) 2 + 8 2 r 12 ] . 
In Fig. (3) we plot the f(r) function versus r depending on different a 3 ^ = a 2 ■ 



(38) 
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VI. THE GENERAL CASE, a 3 / / a 2 



For N > 7 (n > 5), with «3 ^ 7^ a 2 we shall see the role of the third order Lovelock 
parameters as well as the second order. This leads us to a tedious set of differential equations 
which fortunately reduces to Eq. (18) and can be integrated exactly. Indeed, the general 
solution of the Eq. (18) with arbitrary values of a 2 and 0:3, in any arbitrary dimension 
N > 7 can be expressed by the following expression 



/(r) = l + 



(4fi 2 /A)5 



6 (n + 1) (n — 3) a 3 r™ 3a 3 
r n+4 {n - 3) (n + 1) (a| - 3a 3 ) / 16 



6«3 



in which we use the abbreviations 



(39) 



fi = _ r 2n+2 ( n _ 3 )2 (1+n)2>< 

{-3^3 + [-a 2 (2a 2 - 9a 3 ) (1 + n) (n - 3) r 4 + 27 (n 2 - l) r 
+ 27a 2 3 mr 3 (1 + n) (n - 3)} , 



(40) 



5 = 54 (1 + n) z <j ~-a\ (n - l) 2 r 2n + 



[n - 1) (n - 3) a 2 (^3 - -a 2 J r2n+4 + Sa^m (n - 1) (n - 3) r n+3 
+ I ( &3 ~ \ & 2 ) (n _ 3)2 r8+2 " + (* 3 " IN) (n _ 3)2 



+ -alrnr 6 (n - 3) z } , 



(41) 



A = 2 {n - 3) 2 ||v^a 3 + (1 + n) x 



27 _ 9 

— — (n — 1) a 2 — -a 2 a 3 r 4 (n — 3) + c^r 4 (n — 3) ) r T ' 



-—almr 3 (n — 3) 



r 



2n+2 



1 + n) 2 . 



(42) 



To proceed further with these expressions doesn't seem feasible from technical points, 
therefore in the sequel we shall adopt a special simplifying relation between a 2 and a 3 . 
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A. The EYML solution in 7— dimensions with a% = a 2 /3 

In 7— dimensional space time we can see the roles of both second and third order Lovelock 
parameters simultaneously. One of the simplifying, yet interesting case in the solution (43) 

~ 2 

can be obtained if one sets a 3 = The metric function /(r) reads rather simple 



/(r) = l + 



1 - 



^ 3a 2 m 3 (n — 1) a 2 



(43) 



n 2 ] r n+1 ' (n-3)r 4 

which yields an asymptotically flat metric. For r — ■> oo in order to have black hole solutions 
one should investigate the existence of roots of the metric function (i.e., f(r) = 0). To this 
end, one finds the solutions of the following equation 

n — 1 \ n 1 



r + [ a 2 - 



n-3j 3 2 



mr° = 0, 



(44) 



which generally seems a difficult task. But in seven dimensions we get the following roots 



in which for 1 > \/^r an d 



1 - 




(45) 



~ 2 2 ~ 2 

^ — (l — < m < ^y, so that we have both inner and 



outer radii of black hole. In the upper and lower limits of m, we will get just the radius of 
the event horizon of the black hole solutions, i.e., when m — then r + = ^2 — 5t2 while 



^ while r_ does not 



r_ = 0. On the other hand when m — ^ — (l — ^) , r + = yl 
exist. For the choice a 2 = 2, we can define a critical IIlclSS, merit 

= 4/3, so that for m > m crit 

we can have a black hole solution. 

The surface gravity in iV = 7 has the form 



K 



in which 



a 2 



A = 1 + 



1 - 



1 



A( 2 / 3 ) 
3a 2 m 6a 2 



A( 2 / 3 ) 



(46) 



,4 ' 



where r + is the radius of event horizon of the possible black hole. For instance when a 2 = 2, 
(m — l) 1 ^ 4 ) one gets 

\2 



(i.e., r + 



K 



r (3m - 4) 2 [(3m - 4) S 2 / 3 - (3m + 4)] 

9v 3 — -r. 

2 (9m- 12) 11/4 S 2 / 3 



(47) 
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in which 



„ 162m 36 

(9m - 12) 3/2 (3m - 4) ' 



The associated Hawking temperature T# = ^ can be found by using the above result, 
on which one may expect that, T H is strongly a 2 (and <5 3 ) dependent. 



B. Asymptotically dS (AdS) property 



The general solution (39) as r — > oo reads 



/oo M =1 + 



3« 3 y 



1/3 



+ 



3 eta 



1/3 



a 2 
3q;q 



(48) 



where 



3as\j3 (4a 3 — <5f) + 2«2 _ 9a 2 Q ; 3 



It is observed that the metric function (39), can be rewritten as a N (n = N — 2) 
dimensional dS space in which (r) =1 — ( TO „ 3 2 )|„„ 4 ) r2 one defines a cosmological constant 
without cosmological constant— A as an effective cosmological constant given by 



A = 



(n - 3) (n - 4) 



a 



3«3 



1/3 



+ 



1 



3(5.3 



1/3 



«2 
3«3 



(49) 



It is seen that in this effective cosmological constant both a 2 , and 0:3 play role. One can 
easily show that, depending on a 2 and 0J3, A can take zero, positive or negative values, and 
consequently the general solution becomes asymptotically flat, dS or AdS, respectively. For 
instance, from Eq. (43), one can show that the choice <5 3 = o|/3 is asymptotically flat. 

These results verify that, the Lovelock parameters a 2 and a 3 significantly modify the 
properties of EYM black holes as well as their asymptotic behaviors. 



VII. CONCLUSION 



We introduced YM fields through the Wu-Yang Ansatz into the third order Lovelock grav- 
ity with spherical symmetry. The Ansatz and symmetry aided in overcoming the technical 
difficulty and obtaining exact solutions in higher dimensions. In this sense our work is partly 
an extension of our previous work which included only the GB parameter and for N = 5 
[5]. Our solutions include black hole possessing parameters of mass, magnetic charge (which 
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is scaled to Q=l), a 2 (for N > 5) and a 3 (for N > 7). Our analysis indicates that higher 
order Lovelock parameters have less significant contributions compared with the lower order 
terms superposed to the EH Lagrangian. At least this is the picture that we expect in case 
that «fc (k > 4) terms are taken into account. Depending on the choice/relative magnitudes 
of the parameters, formation of the black holes with inner/outer horizons is conditional. 
Asymptotic solutions give rise to dS (AdS) or flat space times which are topologically trivial 
but from the field theoretic point of view they are rather important. 



[1] D. Lovelock, J. Math. Phys. (N.Y.) 12 (1971) 498. 

N. Deruelle and L. Farina-Busto, Phys. Rev. D 41 (1990) 3696. 

G. A. MenaMarugan, ibid. 46 (1992) 4320; 46 (1992) 4340. 

D. G. Boulware and S. Deser, Phys. Rev. Lett. 55 (1985) 2656. 

B. Zwiebach, Phys. Lett. B 156 (1985)315. 

R. G. Cai and N. Ohta, Phys. Rev. D 74 (2006) 064001. 

A. Chakrabarti and D. H. Tchrakian, Phys. Rev. D 65 (2001) 024029. 
[2] M. H. Dehghani and M. Shamirzaie, Phys. Rev. D 72 (2005) 124015. 

[3] T. T. Wu and C. N. Yang, in Properties of Matter Under Unusual conditions, edited by H. Mark 

and S. Fernbach (Interscience, New York, 1969), p. 349. 
[4] P. B. Yasskin, Phys. Rev. D 12 (1975) 2212. 

[5] S. H. Mazharimousavi and M. Halilsoy, Phys. Rev. D 76 (2007) 087501. 

Y. Brihaye, A. Chakrabarti and D. H. Tchrakian, Class. Quant. Grav. 20 (2003) 2765. 
[6] C. W. Misner, K. S. Thorne, and J. A. Wheeler, Gravitation (Freemann, San Fransisco, 1973). 
[7] S. H. Mazharimousavi and M. Halilsoy, Phys. Lett. B 659 (2008) 471. 
[8] S. A. Ridgway and E. J. Weinberg, Phys. Rev. D 52 (1995) 3440. 
[9] E. E. Donets and D. V. Gal'tsov, Phys. Lett. B 352 (1995) 261. 

T. Torii, H. Yajima and K. Maeda, Phys. Rev. D 55 (1997) 739. 



13 



VIII. FIGURE CAPTIONS: 

Fig. (1): The horizon radius r\ versus a.2 for 0:3 = , and different values for mass (each 
mass is written on the correspondence curve) in N=6 and N=7. 

Fig. (2): Plot of the metric function / (r) versus r for N=7, m=l, a 3 — 0, a 2 — 5, 1, 0.5, 
0.1, 0.01 and 0.001. 

Fig. (3): Plot of the metric function /(r) versus r for N=7, m=l, a 2 = 0, a 3 = 5, 1, 0.5, 
0.1, 0.01 and 0.001. 
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